Mem. S.A.It. Vol. 75, 282 
©SAIt 2008 



Radio relics in the MareNostrum Universe 

S. E. Nuza 1 , M. Hoeft 2 , S. Gottlober 1 , R. J. van Weeren 3 , and G. Yepes 4 

1 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 
e-mail: snuza@aip.de 

2 Thiiringer Landessterawarte, Sternwarte 5, 07778 Tautenburg, Germany 

3 Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, Netherlands 

4 Grupo de Astroffsica, Universidad Autonoma de Madrid, Cantoblanco, 28039 Madrid, 
Spain 



Abstract. We identify shocked gas in simulated gala xy clusters extracted from the 
MareNostrum Universe simulation dGottlober et al.l l2006) assuming that shock waves are 
regions of electron acceleration. We perform fl ux number counts within the framework of 
the non-thermal emission model developed bv lHoeft et alj j2008h . Results are presented at 
two different observing frequencies, i.e. 1.4 GHz and 120 MHz, posing interesting con- 
straints for LOFAR and upcoming radio telescopes. 
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1. Introduction 

Radio relics are elongated structures located 
in the outskirts of galaxy clusters. Seemingly, 
these radio features are produced during clus- 
ter merger events as a result of the formation 
of shock waves, which are believed to be ideal 
sites for electron acceleration. In the presence 
of magnetic fields accelerated electron popu- 
lations are capable of producing the observed 
synchroton emission. 

During the last years, the number of 
observed relics has dramatically increased 
due to the spectacular improvement in 
the instruments sensitivity. In particular, 
the system found in the galaxy cluster 
CIZA J2242.8+5301 is one of the most 
spectacular relics known to date, giving 
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strong support to the prev iously mentioned 
as trop hy sic al scenario (Ivan Weeren et al.l 
l201(l . 

In this work, we estimate the amount of 
diffuse radio emission produced in simulated 
relics using a synthetic galaxy cluster sample 
extracted from a cosmological simulation box. 
After identifying shocks in the simulation we 
apply a radio emission model based on the dif- 
fusive shock acceleration (DSA) process and 
count how many relics are observable as a 
function of redshift to assess expectations for 
the upcoming LOFAR radio surveys. 

In Sect. [2] we give a brief overview of the 
method used and in Sect. [3] we close the contri- 
bution with the discussion and comment about 
future work. 
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2. Method 

We select the 500 most massive clus- 
ters i n the MareNo s trum Universe simula- 
tion dGottlober et al.l 120061) . a non-radiative 
smoothed particle hydrodynamics cosmologi- 
cal run with 2 x 1024 3 gas and dark matter par- 
ticles, at 3 different cosmic times (z = 0, 0.5 
and 1). 

Shock waves in the simulated clusters are 
identified in order to compute the non-thermal 
emission as a function of Mach number, mag- 
netic field and thermodynamical state of the 
post-shock region. We assume that the mag- 
netic field in the shocks scales with gas den- 
sity simply obeying flux conservation, i.e. B = 
fioX(n e /«o) 2 ' 3 . where n e is the electron density 
and «o and Bo are two reference values chosen 
to be 10~ 4 crrT 3 and 0. 1 /iG respectively. In our 
model electrons are accelerated by means of 
the DSA mechanism whereas they cool down 
due to synchroton looses and inverse compton 
scattering with the cosmic microwave back- 
ground. For more details about the shock find- 
ing technique and the radio model used see 
lHoeftetalJd2008h . 

For each simulated cluster we estimate the 
radio emission produced in shocks inside a 
sphere of size ~ 2 x R v i r with R v i r being the 
virial radius. In order to estimate the total 
non-thermal emission for less massive haloes 
at each redshift we construct the radio power 
distribution for relics in clusters with differ- 
ent masses and extrapolate the obtained trend. 
From this information we evaluate the proba- 
bility of finding radio relics for a given clus- 
ter mass, radio power and cosmic time. To 
compute the radio relic luminosity function 
(RRLF) as a function of redshift we convolve 
this proba bility with the halo m ass function 
given bv ISheth & Tormenl (120021) . Finally, us- 
ing the RRLF it is possible to perform the flux 
number counts. 



3. Discussion 

Figure [TJ show the model predictions for the 
observed relic cumulative number per logarith- 
mic radio flux at 1.4 GHz (normalized to the 
NVSS sample at 100 mJy) and 120 MHz (cor- 
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Fig. 1. Cumulative number of observed relics as a 
function of radio flux. The histogram shows the ob- 
served NVSS relic sample (S > -40°) as a reference. 
Also shown are models at 1.4 GHz (NVSS) and 120 
MHz (LOFAR Tier 1). 



responding to the LOFAR Tier 1 survey). It can 
be seen that LOFAR will be able to observe 
thousands of relics. The actual number will de- 
pend on the final sensitivity of the survey and 
on the details of the shock magnetic field in- 
volved. We will present a comprehensive anal- 
ysis using different model assumptions in a 
forthcoming paper (Nuza et al. 201 1, in prepa- 
ration). 
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